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ABSTRACT 

This  report  is  concerned  with  the  development  of  a  new  approach 
to  the  problem  of  stability  for  multidimensional,  causal,  recursive, 

*all  pole*,  digital  filters.  The  distinguishing  feature  of  this  approach 
is  that  general  stability  criteria  can  be  derived  directly  in  terms  of 
the  coefficients  of  the  transfer  function  of  the  filter.  Thus  by  use  of 
this  method  it  is  sometimes  possible  to  determine  which  coefficients 
of  the  transfer  function  are  critical  to  the  stability  of  the  filter, 
information  which  is,  of  course,  important  in  filter  design.  Also  the 
emphasis  of  this  approach  is  on  the  development  of  a  conceptual  method 
for  considering  the  problem  in  complete  generality. 
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The  Coefficient  Space  Approach  to  the  Stability 
of  Multidimensional  Digital  Filters 


1.  Introduction 

Consider  a  filter  whose  transfer  function  has  the  form 

H(z)  =  1/(1  -  Q(z‘b)  (1.1) 

where  z  =  z^)  is  a  vector  complex  variable,  and  QCz)  is  a 

polynomial  in  N  variables  with  real  coefficients  and  zero  constant  term. 
Such  a  transfer  function  describes  an  ’all-pole’  ’causal’,  recursive, 
multidimensional  digital  filter,  and  it  is  known  that  questions  of  stability 
generally  reduce  to  a  question  of  stability  for  filters  of  this  type.  It 
is  also  known  that  such  a  filter  is  stable  iff 

P(z)  =  1  -  Q(z)  (1.2) 

has  no  zeros  S  =  •••>  S^)  such  that  |S^|  ^1  for  all  i  .  For 

this  reason,  it  is  convenient  terminology  to  state  that  a  polynomial  P(z) 
of  several  variables  is  stable  if  P(S)  =  0  implies  that  |b^|  >  1  for 
some  i  ,  The  question  of  stability  then  becomes,  primarily,  the  problem 
of  determining  if  such  a  polynomial  is  stable. 

An  N-tuple  of  the  form  jf  =  (j^,  ...,  where  each  is  a  non¬ 

negative  integer  is  called  a  vector  index.  If  for  such  an  index  X  > 
define  ‘z’^  by 


then  every  polynomial  P(z)  in  N-variables  can  be  written  in  the  form 


1 


(1.4) 


P(z)  =  ^  ^ 

3 

where  for  each  ^  ,  is  a  real  numberj  equal  to  zero  except  for  a 

finite  number  of  indices  (see  Davis  and  Souchon  [1975]).  For  example  if 

2 

=  1.00  +  O.Sz^z^  -  .llz^z^ 


then  a^Q  =  1.00,  a^^^^  =  .05,  =  -.11  and  all  other  coefficients  equal 

zero. 

Suppose  that  a  fixed,  ordered  set  (j*. ,  ...,  "j"  )  of  non-zero, 

jL  Z  K 

N-dimensional  vector  indices  is  given,  and  consider  the  set  of  polynomials 
of  the  fomn 


P(^  =  1 


(1.5) 


where,  for  each  i  ,  a^  is  a  real  number.  Each  such  polynomial  (1.5) 

■^i 

is  naturally  associated  to  the  point  (aj.  ,  . . . ,  a^  )  of  K-dimensional 

^1 

euclidean  space.  Conversely,  each  point  (A^^,  •••, 
euclidean  space  is  associated  to  a  polynomial  of  type  (1.5)  by  the 
correspondence 


A^)  of  K-dimensional 


A «  ~  a  ^  j  x~l,..«,K  (1*6) 

relative  to  the  ordered  set  of  indices  .  •  •  ,  •  Thus  relative  to 

i  K 

this  ordered  set  of  indices,  K-dimensional  space  becomes  the  coefficient 
space  for  polynomials  of  type  (1.5).  A  point  (A^ 

coefficient  space  will  be  said  to  be  a  stable  point  if  its  associated 
p  o lynomial ,  namely 


of  this 


2 


(1.7) 


P(t)  =  1  -  XI  ^ 
i=l 

is  stable.  The  set  of  all  such  stable  points  in  the  coefficient  space  will 
be  called  the  region  of  stability.  The  problem  that  will  be  considered 
here  is  what  can  be  found  concerning  the  region  of  stability  for  a  given 
set  of  vector  indices,  or  equivalently,  for  a  given  type  of  polynomial. 

In  certain  cases,  the  region  of  stability  can  be  specified  completely. 
Consider  the  ordered  set  of  indices  =  (10),  'j*'2  =  (01)  “  (H)*  The 

associated  polynomials  are  of  the  form 

P(z^,Z2)  =  1  -  -  k^z^  -  k^z^z^  ,  (1.8) 

and  the  associated  coefficient  space  is  3  dimensional  euclidean  space. 

Huang  [1972]  has  shown  that  (1.8)  is  stable  iff 

IA3I  <  1 

IA1+A2I  <  1-  A3 

|a^  -  A2I  <  1  +  A3  . 

The  region  of  stability  is  illustrated  in  Figure  1.  Also,  it  can  be 
shown,  see  Jury  [1974]  that  a  similar  type  of  specification  can  be  given 
for  any  region  of  stability  associated  to  polynomials  in  one  variable,  or 
equivalently,  to  any  ordered  set  of  1-dimensional  vector  indices. 

Although,  regions  of  stability  can  be  very  difficult  to  determine  for 
more  complicated  types  of  multivariable  polynomials,  in  the  following  it 
will  be  shown  that  several  properties  of  regions  of  stability  can  be 
derived  in  general. 


3 


As 


FIGURE  1 
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2. 


The  Basic  Theorems 


Throughout  the  following  it  will  be  assumed  that  "J. ,  . . . ,  ^  is  a 

X  K 

given,  fixed,  ordered  set  of  N  dimensional  vector  Indices.  Thus  the 
correspondence  between  polynomials  of  the  form  (1.7)  and  points 
(Aj^,  . . . ,  Aj,)  of  K-dimensional  space  is  assumed  fixed. 

The  order  of  a  vector  index  "f  =  (j^,  ...,  j^)  is  defined  by 

order  j  =  j  ^^  +  . . .  +  . 

A  coefficient  A^  of  a  polynomial  of  type  (1.7)  is  called  a  leading 
coefficient  if  its  associated  vector  index  "t.  has  maximal  order.  A 

-  -’i 

polynomial  of  several  variables  may  have  several  leading  coefficients. 

Theorem  2.1  If  (A^,  . .  . ,  A^,)  is  a  stable  point  of  the  coefficient  space 
and  if  A  is  the  sum  of  the  leading  coefficients  of  the  polynomial 
corresponding  to  this  point,  then  |a|  <  1. 

Proof.  Let  PC?)  be  the  stable  polynomial  associated  to  (A^,  ..., 

It  follows  that  the  one  variable  polynomial  p(z)  defined  by 

p(z)  =  P(z,  ...,  z)  (2.1) 

is  also  stable.  Moreover  the  leading  coefficient  of  p(z)  is  seen  to  be 
A  .  The  polynomial  p(z)  factors, 

p(z)  =  A(z-t^) . . . (z-t^)  (2.2) 

where  t^,  ...,  t^^j  are  its  complex  roots.  Since  p(z)  is  stable  it 
follows  that 

|t^|  >1  i  =  1,  ...,  M  .  (2.3) 
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But  then. 


1  =  |p(o)|  =  |a| 


M 

I  n 


i=l 


(2.4) 


where 


M 

I  n  t  I  >  1  . 
i=l  ^ 

Hence, 


1a|  <  1  . 

Theorem  2.2  (Necessity  Theorem)  Let  ct  =  (a^^, 
of  complex  numbers  such  that  for  each  i  =  1, 
for  each  i 


s 


i 


(2.5) 


(2.6) 


,  Oj^)  be  any  N-tuple 
N  ,  I I  =  1  ,  and 


(2.7) 


is  a  real  number.  Then,  if  (A^^,  . . . ,  A^)  is  any  stable  point,  it  must 
be  true  that 


SiAi+  ...  +s^<  1  .  ,  (2.8) 

(Note  that  since  each  s^  has  modulus  1  and  is  real,  each  s^  must  equal 
+1  or  -1.) 

Proof.  Let  PC^  be  the  stable  polynomial  associated  to  (A^^,  ..., 

Suppose  a  satisfies  the  hypotheses  and  define  p(z)  by 


p(z)  =  ^(“iZ,  . . . ,  a^z)  . 


(2.9) 


Then  p(z)  is  a  one  variable  polynomial  with  real  coefficients.  Moreover 


6 


p(l) 


=  1  - 


K 

z 

i=l 


A^(a^, 


K 

A. s ,  . 
i  1 

i=l 


(2.10) 


If 


K 


E 


SiAi  1  1 


(2.11) 


then 


p(l)  1  0  . 


(2.12) 


But  p(0)  =  1,  and  p(z),  as  a  function  of  a  real  variable,  is  clearly 
continuous.  Therefore  by  the  intermediate  value  theorem  for  real  functions, 
there  must  exist  a  real  number  t  ,  0  <  t  1  ,  such  that  p(t)  =  0.  But 
then 

P(aj^t,  ...,  =  0  (2.13) 

and  for  each  i  =  1,  ...,  N  ,  ^  1  ,  a  contradiction  to  the  assumption 

that  PCz)  is  stable.  Therefore  it  must  be  true  that 


.A. 

1  1 


< 


1  . 


(2.14) 


Theorem  2.3  (S3niimetry  Theorem)  Let  a  =  (a^,  .  . .  ,  a^)  be  any  N-tuple  of 
complex  numbers  such  that  for  each  i  =  1,  ...,  N  ,  |a^|  =  1  ,  and 


s . 

1 


s  a 


(2.15) 
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is  a  real  number.  Then  a  point  (A^,  . . . , 

stable.  (Note  as  before  that  each  equals  +1  or 

-1.) 

Proof.  Let  P(Zj^,  ...,  be  the  polynomial  associated  to  (A^^,  ...,  • 

Define  P*  (z^^,  ...»  z^^)  by 

P*(Zj^,  ...)  z^)  =  P(aj^z^,  •.«>  *  (2.16) 

It  is  not  difficult  to  check  that  P'  is  the  polynomial  associated  to 

P'(6j^,  ....  =  0  iff  P(8j^j  ...»  8^)  =  0  where  8j^  =  and  for 

each  i  |6!i  =  |a^8^|  =  since  ia^|  =1  . 


. ..,  Sj^Aj,) .  Moreover  P'  is  stable  iff  P  is.  For 


A^)  is  stable  iff  the  point 


Theorem  2.4  (Sufficiency  Theorem)  Let  (A^^,  . . . ,  A^)  be  a  point  such 
that 


(2.17) 


Then  (Aj^^  . . .  ^ 

Proof.  Consider  the  polynomial  (1.7) 


A^)  is  stable. 


P(z^, 


If  P(8)  =  0  where  |8^|  ^1  >  then 


(2.18) 


(2.19) 


Therefore 


1  =  E  A.  <  E  l^il  <  1 

i=l  i=l 


(2.20) 


a  contradiction. 
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Corollary  2.1  Let  . A^,)  be  a  point  such  that  A^  ^  0  for  all 

X  •  Then  ^A^^j  • « •  > 

K 

X)  \  1  •  (2.21) 

i=l 

Proof.  Apply  Theorem  2.2  and  Theorem  2.4. 

Theorem  2.4  and  its  corollary  can  be  given  simple  geometric  inter¬ 
pretations.  The  region  of  points  (A^,  ..., 

1Ai1  +  ...  +  |a^1  <  1  (2.22) 

describe  the  K-dimensional  'diamond',  centered  at  the  origin  and  whose 
points  lie  along  the  axes.  Theorem  2.4  states  that  the  K-dimensional 
diamond  is  wholly  enclosed  by  the  region  of  stability.  The  Corollary 
states  that  in  the  positive  'quadrant',  the  region  of  stability  always 
coincides  with  the  'diamond'  (see  Figure  1).  In  the  following,  the 
Symmetry  Theorem  will  be  applied  to  show  that  regions  of  stability  also 
satisfy  certain  geometric  symmetries. 

3.  Symmetry 

As  the  proof  of  Theorem  2.3  shows,  symmetries,  (in  this  case  sequences 
of  reflections)  between  the  points  of  the  region  of  stability  arise  from 
transformations  of  the  variables  of  the  associated  polynomials.  The 
transformation  of  variable  is  given  by 

(z^,  . . . ,  •  *  •  >  ~  . . . ,  cijjZjj)  (3. 1) 

where  (a^,  ...,  a^)  is  a  complex  vector  satisfying  the  conditions  that 
for  each  i  ,  |a^|  =  1  and 


A^)  is  stable  iff 
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a 


(3.2) 


s .  = 


is  a  real  nuinber.  Such  a  stability  invariant  transformation  transforms  the 
coefficients  (A^, 

always  equals  +1  or  -1,  this  transformation  is  geometrically  interpreted  as 

a  sequence  of  reflections  of  the  axes. 

N 

At  least  2  vectors  a  can  be  obtained  by  choosing  each  equal 

to  +1  or  -1.  The  symmetries  of  the  coefficient  space  thus  obtained, 
however,  are  not  necessarily  distinct,  and  may  not  include  symmetries  which 
can  be  obtained  by  allowing  the  to  be  complex.  However,  symmetries 

so  obtained  are  more  easily  studied  and  for  this  reason  Xi;ill  be  called 
simp le  symmetries.  In  the  following  we  will  restrict  our  attention  to  the 
study  of  simple  symmetries. 

The  transformation  (symmetry)  of  the  coefficient  space  determined  by 

the  vector  "a  can  be  described  completely  by  the  vector  s  =  (s^,  ...,  s^) 

X  K 

specified  by  eqn.  (3.2).  Each  coordinate  of  the  vectors  "a  =  (a^,  ...,  ot^) 
and  ^  =  (s- ,  ...,  s^)  equals  +1  or  -1.  However  for  reasons  which  will 
become  clear,  it  will  be  more  convenient  to  use  1  in  place  of  -1,  and  0 
in  place  of  1.  The  multiplicative  relation  (eqn.  3.2)  between  the  vectors 
's  and  ^  now  becomes  the  following  additive  one  in  modulo  2  arithmetic. 


into 


Since  each  s . 


=  "a  •  (mod  2)  (3.3) 

where  the  * • '  represents  the  vector  dot  product  of  the  0,  1  vector  "a 
with  the  integer  vector  index  .  If  the  N  x  K  integer  matrix  J  is 
defined  by 

J  =  (t/.  ....  ij)  (3.4) 
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then  the  vector  s  which  describes  the  simple  symmetry  arising  from  the 
change  of  variables  described  by  the  vector  a  satisfies 

^  J  (mod  2).  (3.5) 

N 

Moreover  the  2  possible  choices  for  a  can  now  be  viewed  as  the 
elements  of  the  N-dimensional  vector  space  over  the  Galois  field,  GF(2), 
of  two  elements,  a  tool  familiar  in  algebraic  coding  theory  (Berlekamp 
[1968]).  The  simple  symmetries  can  now  be  easily  classified  using  the 
linear  algebra  of  these  vector  spaces. 


Theorem  3.1  The  set  of  K-dimensional  vectors  over  GF(2)  which  describe 
the  set  of  simple  symmetries  of  the  region  of  stability  is  the  set  of 
vectors  spanned,  modulo  2,  by  the  row  vectors  of  the  matrix  J  . 

Proof.  Each  vector  "s  arises  by  eqn.  3.5  from  an  "a  .  Each 
a  =  (a^,  ...,  Oj^)  can  be  written  as 

a  =  a^e'j^  +  ...  +  a^^-e^  (3.6) 

where 

SjL  ~  C0>  •  •  •  >  1>  •  •  •  >  0)  (3*  7) 

with  a  1  in  the  i^^  coordinate.  But  then 


■^  =  a^(^j^J)  +  ...  +  a^^Ce^^J) 


by  linearity;  and  for  each  i  ,  "e^J  is  the  i^^  row  of  J  . 


(Note  also  that  by  the  above  theorem  the  set  of  simple  symmetries  forms 
an  abelian  group.) 
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Corollary  3.1  The  number  of  simple  symmetries  equals  2^  ,  where  L  is 
the  modulo  2  rank  of  the  matrix  J  . 

Proof.  Immediate. 

Example  3.1 

Consider  the  filter  (eqn.  1.8)  studied  by  Huang.  The  vector  indices 
in  this  case  are 

=  (10) 

^2  =  (01)  (3.8) 

=  (11)  . 

The  matrix  J  is  therefore 


J 


1  0 
0  1 


(3.9) 


2 

The  modulo  2  rank  of  J  is  clearly  2.  Thus  there  are  4  (=2  )  simple 
symmetries.  Each  symmetry  is  described  by  an  element  of  the  row  space  of 
J  .  They  are: 


=  (000) 

S2  =  (101) 

i*3  =  (011) 
i*4  =  (110) 

which  describes  the  symmetries 

(A^jA2>A2)  (A^jA2>A2) 

(~Aj^ J  A2  >  “A^) 
(Aj^ )  *A2  >  “A^) 
(~Aj^>~A2> A^)  * 


(3.10) 


(3.11) 
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Note  that  in  Figure  1,  there  are  only  two  different  basic  shapes  for  the 
region  of  stability  in  the  eight  different  quadrants  of  the  coefficient 
space,  and  that  consistent  with  the  above  symmetries  each  shape  occurs 
symmetrically  in  four  quadrants. 

Example  3.2 

Consider  a  filter  whose  transfer  function  is 


(3.12) 


where 


-  k^z^z\z'^^z^  (3.13) 


In  this  case  the  vector  indices  are 


=  (1301) 
(1221) 


(1331) 


(3.14) 


3^4  =  (0331) 
T5  =  (1031) 


And  the  matrix  J  is  therefore 


J 


0  2  3  3  : 


1110] 


3  2  3  3  C 


1111] 


(3.15) 


Calculating  modulo  2,  and  row  reducing 
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J 


(3.16) 


1 

0 

0 

1 


1 

1 

1 

1 


0 

1 

1 

1 


-> 


1  1 
1  0 
0  1 
0  0 


0 

1 

1 

1 


Therefore  J  has  rank  4,  and  there  are  2  =  16  possible  distinct  simple 

symmetries.  Moreover  each  possible  symmetry  can  be  described  by  a  modulo 
2  sum  of  the  rows  of  J  .  For  example  adding  every  row  we  obtain 


■s  =  (1  0  0  1  1) 


(3.17) 


which  corresponds  to  the  s3nnmetry 


(Aj^» A2> Ag)  (“A^> A2 » A^j ~A^ , “A^)  (3.18) 

of  the  coefficient  space.  Thus,  for  example,  since  we  know  that  in  the 
positive  quadrant  the  shape  of  the  region  of  stability  is  the  part  of  the 
diamond  in  that  quadrant  it  follows  that  in  the  quadrant  which  corresponds 
to  the  above  symmetry,  the  region  of  stability  is  again  the  part  of  the 
diamond  in  that  quadrant. 

Example  3.3 

Consider  the  class  of  polynomials  above  without  the  last  term. 

=  1  -  Aj^z^Z2  -  k^z^z^z^z^  -  k^z^z^z^z^ 

-  A^z^z^z^  .  (3.19) 

Since  the  rank  of  J  will  still  be  four,  and  the  dimension  of  the 
coefficient  space  is  four,  it  follows  that  every  quadrant  is  symmetrical 
to  every  other  by  an  appropriate  simple  symmetry.  That  is,  every  possible 
change  of  sign  will  occur  among  the  simple  symmetries.  It  follows  that 
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the  region  of  stability  is  the  diamond,  and 


<  1 


is  a  necessary  and  sufficient  condition  for  the  stability  of  polynomials 
of  this  type. 
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